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I. INTRODUCTION 



The experimental realization of Bose-Einstein conden- 
sation in atomic or molecular gases has rekindled interest 
in the quantum theory of dilute bosonic gases JlJ, |2, 13, IJ] . 
So far, most investigations were based on the s-wave scat- 
tering approximation which neglects microscopic details 
of the two-particle interaction (such as attractive regions 
of the potential), and replaces the true interaction by 
a contact interaction pseudopotential. To describe the 
impact of a possibly complicated microscopic structure 
of the potential on quantities like the occupation num- 
ber of the condensate, it is however necessary to go be- 
yond the s-wave approximation. On the other hand, the 
current theoretical description of Bose-Einstein conden- 
sates is predominantly based on a mean-field expansion 
around a macroscopically occupied state representing the 
condensate. In leading order this expansion yields the 
Gross-Pitaevskii equation, which forms the basis of most 
theoretical approaches to Bose-Einstein condensed gases 
i 

The scope of the present article is to investigate con- 
sequences of deviations from pure s-wave scattering for 
the applicability of this mean-field expansion, and to de- 
rive the general zero-temperature equations governing 
the condensate and quantum fluctuations above it. To 
this end, we present an ah initio approach, controlled 
by properties of the bare interaction potential in the mi- 
croscopic second-quantized equation of motion. While 
the presented approach is essentially restricted to suffi- 
ciently "weak" , that is, integrable potentials, we do not 
use pseudopotential or T-matrix approaches 0, Q, 
for which the direct connection to microscopic proper- 
ties of the interaction potential is lost. Our approach 
can thus be viewed as being complementary to the latter 
treatments. In particular, if one can separate the po- 
tential into a strong ultra-short-range part, assumed to 
be replaceable by an integrable (pseudo-)potential, and 
a weaker remaining contribution with longer range, our 
analysis applies to this latter part as a low-energy, effec- 
tive description. 



We test the validity of the mean-field expansion, which 
is usually an ad hoc assumption, employing well-defined 
scaling properties of the interaction potential with the 
number of particles and thus the density of the system, 
under the prescription that the system always remains di- 
lute. It turns out that the consistency of the mean-field 
expansion sensitively depends on the scaling exponent 
of the interaction potential, and thus on the (formal) de- 
pendence of the interaction range on the particle number. 
This formal dependence on the particle number can be 
made more intuitive by comparing the dilute-gas limit 
with the thermodynamic limit. 

The structure of the paper is as follows. In the next 
Section we derive, using a number-conserving mean-field 
expansion, the equations of motion for the mean field it- 
self (a modified Gross-Pitaevskii equation) and for the 
single- and multi-particle excitations above mean-field. 
In the third Section, the consistency of this mean-field 
expansion is tested with three scaling behaviors of the 
particle interaction potential. The fourth Section treats, 
as an explicit example, a dipole-dipole interaction poten- 
tial. Even though such a potential is on the borderline 
of applicability of our approach [because it is (just) not 
integrable in three spatial dimensions] , we show that the 
calculation of certain quantities like the quantum deple- 
tion of the condensate is still feasible. 



II. NUMBER-CONSERVING MEAN-FIELD 
EXPANSION 

We assume in the following that the size of the inter- 
acting atoms or molecules forming the dilute gas is much 
smaller than all other involved length scales such as the 
s-wave scattering length a^; typical energy scales asso- 
ciated with scattering processes are thus assumed to be 
well below the internal energy differences. In this limit, 
we may adopt the point-particle approximation and the 
interacting objects can be treated as structureless enti- 
ties described by the field operator ^. Setting Planck's 
constant h and the mass m of the atoms or molecules 
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equal to unity, the Heisenberg equation of motion for the 
field operator reads 

^^ - -- + K.(r)U(r) 



dPr' i>\r')Vint{r - r')*(r')*(r-) , (1) 



with Vcyit{i') denoting the external one-particle poten- 
tial of the trap. The finite-range two-particle interac- 
tion potential 14nt('" — is assumed to be sufficiently 
"weak" in the sense that the Born approximation ap- 
plies, e.g., Vi-nt{r) is supposed to be integrable (i.e., to 
decay more rapidly with distance than such that its 
Fourier transform 14nt(fe) is a well-defined function of k. 
For large wavelengths, the finite-range two-particle inter- 
action potential Viat{r — r') yields an effective coupfing 
constant g determined by the s-wave scattering length ag 
(in three spatial dimensions) 



g 



47ra, 



dVT4„t(r) = T4„t(fc-0) 



(2) 



The Bogoliubov mean-field approximation U is based on 
the assumption that the fluctuations Xj describing single- 
particle excitations above the ground state are small 
(compared to the mean field) if the number N of parti- 
cles is large. In the following, we present a formal deriva- 
tion of the mean-field expansion for the system described 
in Eq. in the large-A^ limit, i.e., an expansion of ^ 
into powers of N . We shall employ a particle-number- 
conservmg mean- field ansatz (satisfying (#) ^ 0) [lllIT^ 



+ X + C 



A 



(3) 



Here, the order parameter is i/'o- The single-particle ex- 
citations are described by x, where single-particle here 
means that the Fourier components of x ^^'^ linear su- 
perpositions of annihilation and creation operators of 
quasiparticles and a\,, cf. Eq. ||HJ). The higher-order 
corrections C are due to multi-particle excitations and 
correlations. Here N — A^A counts the total number of 
particles. 

The dilute-gas limit may formally be defined by plac- 
ing a very large number iV | cx) of identical bosons into 
a finite trap with gN remaining constant [11 m. The 
diluteness parameter of the gas then scales as (g'^gY^^ oc 
1/7V, i.e., the gas becomes rapidly more dilute with in- 
creasing particle number in this limit. This limit should 
be compared and contrasted with the usual thermody- 
namic limit, in which the density and particle interaction 
remains constant, while the size of the (trapped) system 
increases with TV — > oo, adjusting the trapping potential 
correspondingly (in D spatial dimensions, the thermo- 
dynamic limit corresponds to keeping Nu!^ constant for 
N oo, where lu is the geometric mean of the trapping 
frequencies @). In the presently used dilute-gas limit, on 



the other hand, the trapping potential remains constant, 
but the interaction and the density change. The advan- 
tage of the limit gN constant is that in this limit we have 
a well-defined prescription for checking the mean-field ap- 
proximation, keeping one power of g for each factor of N, 
cf. 2]. 

Given the dilute-gas limit, there arises a question 
which is central to the applicability of the mean-field 
expansion: Do the supposedly small multi-particle cor- 
rections actually become small in this limit for a given 
behavior of the interparticle potential Vintir — r')? In 
order to tackle this question, let us consider the fol- 
lowing scenario: Initially, the interaction is completely 
switched off, Vint{f — f') ^ Q, and all N bosons oc- 
cupy the same single-particle state (at zero temper- 
ature), described by the single-particle wave- function 
tJjo/Vn. In this case, Eq. 131 is trivially satisfied with 
V-o = 0{VN), X - OiN°), C = 0, and i = aofio^^^N^^^, 
where fiQ = ajao corresponds to the macroscopically oc- 
cupied state, cf. jlL] . If we now switch on the interaction 
slowly enough, the system stays in its ground state due 
to the adiabatic theorem and we may follow the evolution 
of the field operator 4' and hence ipo, x, and ( (with A 
remaining constant, i.e., the evolution at given particle 
number) via insertion of Eq. © into (Q) . 

Assuming that the corrections C indeed become small 
for t oo, the leading order (\/A) yields, in the dilute- 
gas limit, the Gross-Pitaevskii (CP) equation 9] for the 
order parameter ij^o 



.gl^ol V'o 



(4) 



The next-to- leading order (A^*^) terms govern the evolu- 
tion of the fluctuations x via the nonlocal analogue of 
the Bogoliubov-de Gennes equations ,lt (lOi] 



■ dx 



= I d\'VUr-r')[\M^x{r')+ijlxHr')\ 



(5) 



The remaining equation of motion for ^ takes for a gen- 
eral nonlocal potential the rather complicated form 



dt V ? + ^- + I' 



+ ^nt* 



H-V'oX^^int *X + ^oiVint * X^)X + "00 Vint * (X^X) 



K„t*(x^x) x + o i-^ 



(6) 



where Vint* is an abbreviation for the convolution with 
the quantity right from the "*" as in Eq. (0, and 
©(C/a/A) denote formally sub- leading terms containing 
(■ or C,^ or both. Observe that the commutation re- 
lations [x, AA"i/2] = [C,AA'-i/2] = vahd initially are 
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preserved under the evolution given by the equations 
above, i.e., the excitations \ and C are always particle- 
numb er- conserving, cf. mill- 

If we start with free particles, and switch on the in- 
teraction by following the evolution in Eqs. ©, and 
© , we observe a mixing of x and due to Eq. ^ , and 
C does not vanish anymore due to the four source terms 

X^^nt * X, {Vint * X^)x, Vint * (x^ x) , and [Vint * {x''x)]x- 
The viability of mean-field theory depends on the scal- 
ing of the induced terms in (■ with particle number, i.e., 
if and to which extent they decrease for t oo. Naive 
power counting would imply ( = 0{1/\/N), but the non- 
linearity due to the product of two or more field operators 
(such as in the above source terms) and the associated 
mode sum(s) may compensate the smallness of 1/\/N. 
As an example, let us consider the expectation values of 
the four source terms. The last term \Vint * (x^x)]x is 
both sub-leading and odd so that its expectation value 
vanishes. However, the expectation values of the remain- 
ing terms do not vanish, and affect {({N)). 

For the sake of notational convenience, we con- 
sider in what follows a homogeneous condensate 
V'o = ^/Qo = const in a constant "trapping" potential 
V^xt = ~gQOi so that the two last terms in the round 
brackets on the right-hand side of the Bogoliubov-de 
Gennes equations © cancel. A normal-mode expansion 
for the fluctuation operators then yields in Fourier space 

^ ( T ^ SoVint{k) \ Xk + QoVnAk)x\ , (7) 



This faster-than-polynomial decrease renders all relevant 
expectation values finite. For example, the expectation 
value of the so-called "anomalous" term reads 



(xMnt * X> 



27r2 



dk 



2ujk 



(10) 



where we have assumed a spherically symmetric potential 
Vi-atiT) and hence V-mt = Mnt(fc)- The expectation value is 
taken in the ground state of the quasiparticles (remember 
the adiabatic switching process) which is annihilated by 
Sfc (but not by Xk, of course). For a sufficiently regular 
14nt(fc), an expansion into powers of l/(CA:cut) yields 



(xMnt * X) 



ii 

27r2 



dk V:i 



9^ 



O 



9^ el 



.(11) 



The l/(^A:cut)-corrections depend on the explicit form of 
Vint{r) and hence Vint{k), e.g., for yi nt(fc) = gB (fccut - k), 
the integral in Eq. ^ yields g^{y/£.'^k^^^ + 4 - 2)/^. 

In contrast to the so-called "anomalous" term 
(x^nt * x) J the expectation value of the quantum de- 
pletion terms V^nt * (x^x) and ((V^nt * X^)x) occurring 
as source terms in Eq. ^ do not have a contribu- 
tion linear in kcut- Hence their (for ^fccut ^ 1) dom- 
inant terms are independent of kcut and both give 
<7(X^X) « y/g^/iSn'^). The higher-order l/i^hnt)- 
corrections, however, again depend on the shape of 
Vint{k) and can be calculated analogously. 



so that the annihilation operators Xk of the original 
bosons have the following Bogoliubov transformation 
form in terms of the quasiparticle operators 0^,0,^: 



Xfe = 



2ujk 



fe2 



flfc 



(8) 



The generalized Bogoliubov dispersion relation reads 



fe^6'oMnt(fc) 
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(9) 



For wavenumbers fc^^^ ^ i much larger than the 



inverse of the healing length ^ 



we have 



Xk ^ o,k — d]^goVintik)/k^ , i.e., the quasiparticles be- 
come asymptotically equal to the original bosons. How- 
ever, this 1/k^ decrease alone is not sufficient for render- 
ing the expectation value of (x^nt * x) finite (except in 
one spatial dimension). 

A smooth (and integrable) interaction potential 
Vint{r) implies a faster-than-polynomial decrease 
Vint(fc^ ^ fccut) ^ for large wavenumbers k^ ^ k^^^. 
The cut-off wavenumber fccut, then, is determined by 
the first significant deviations of Vint from its long- 
wavelength behavior 14nt(fe^ ^ ^cut) = 9 and is assumed 
to be much larger than the inverse of the healing length. 



III. CONSISTENCY OF MEAN-FIELD 
EXPANSION 

To explicitly address our principal question of whether 
the higher-order corrections C in the mean-field expan- 
sion Q indeed become small in the large-iV limit (in 
which one would naively expect mean-field to become 
more and more accurate), we now discuss three concrete 
examples for the formal scaling behavior of the interac- 
tion Vint with the number of particles N. To this end, the 
particle-number-dependent potential V^n^'' is always cho- 
sen such that it reproduces the g cx 1/N dilute-gas limit 
prescription for the behavior of the coupling constant in 
the large- limit. 

In the first example, we assume kcut to be (formally) 
independent of N 



(12) 



Since fccut remains constant for N 00, the mode sums 
due to the nonlinear terms in Eq. © cannot compensate 
the smallness of the pre-factors such as 1/ ^/N and thus 
we obtain 



' /]v 



(13) 
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However, this decrease is not sufficient yet for ensuring 
the usual split of the total density 

Q = (^¥^) - + {xH) + o (^-^) , (14) 

since ?/'o(C) could be of the same order as (x^x)- For the 
above split, we need (C) = 0{l/N) instead of 0(l/\/]V), 
which requires absorbing the expectation values of the 
source terms into a modified GP equation 

= (-^ + Kxt + .<7l^ol'+.g(xtx)) V-o 

+ (xM„t * X>V'o + ((^^nt * X^)x)^o • (15) 

Inserting Eq. (|ll|l . we see that the dominant (for 
^.kcnt S> 1) and density-independent fluctuation contri- 
bution (x'^nt * x) proportional to fccut can be absorbed 
by a renormalization of the coupling constant 

5ren = 9 - ^ J dk V^^^ 

= 5 + J d^rVi^tV-^Vint, (16) 

where the second line is due to Parseval's theorem, us- 
ing the formal inverse of the Laplace operator V^. Since 
— V~^/2 and Vint/'2 are the kernels of the Schrodinger 
Hamiltonian Hq and the interaction Hamiltonian Hi , re- 
spectively, this just corresponds to the usual one-loop 
renormalization of the interaction potential. Note that 
in three spatial dimensions, sub-dominant contributions 
like (x^x) cannot be absorbed by such a renormalization 
procedure. 

The physical significance of the rcnormalized coupling 
f^ron can be demonstrated further by calculating the fccut- 
corrections to the total energy of a homogeneous gas 

N N P ~ 
E = {H)^ TVKxt + ^gg ~^ej dk v;^,, , (17) 

which gives E = NVcxt + A^.grcnf?/2 plus higher-order 
corrections (cf. 2]). 

A potential complication which arises from the pre- 
scription H12|l is that the range of the two-particle in- 
teraction potential in Eq. (|12l) will exceed the inter- 
particle distance di for large N since go ^ N and thus 
di ~ 1/A^^/'^ (in three spatial dimensions). If the two- 
particle interaction potential is mainly caused by direct 
collisions, this might lead to a conflict with the point- 
particle approximation used in writing down the starting 
point of the analysis, Eq. Therefore, as a second ex- 
ample, consider the support of Vint(^) to be decreasing 
in proportion to the inter-particle distance di, 

vLVir) = V^^{N'/'r). (18) 

As a consequence, the cut-off scales as fccut ~ N^^'^, and 
thus (x^^nt * x) ~ N~'^/'^ , which ensures that the remain- 
ing C corrections are still small, though decreasing with 



a smaller power in N than in the first example Eq. (|13|l : 

c = o(^)- (IS) 

The above scaling with particle number is due to the fact 
that for each additional operator in Eq. © , which might 
contribute a factor fccut N^^^ after the mode summa- 
tion, there is a pre- factor of order 1 / ^/N. In this case, 
the split in Eq. H14|l is still possible provided the mod- 
ification of the GP equation (|15|) is employed, but the 
estimate of the accuracy is now the rather slow decrease 
of 0(l/iVi/6) instead of 0(l/iVi/2). 

Finally, as a third example, we investigate a scaling 
employed by Lieb et al. \13L [Tj| , used in a proof of the 
asymptotic exactness of the Gross-Pitaevskii energy func- 
tional in three spatial dimensions, and in a analysis of 
one-dimensional systems of bosons in 3D traps, respec- 
tively. This scaling reads 

vL^\r) = N^V^,\Nr), (20) 

and implies that the cut-ofF increases linearly with par- 
ticle number, fccut N, and thus (xMnt * x) so 
that the anomalous term becomes of the same order as 
other terms in the GP equation, for example as large as 
the mean- field interaction term g|?/'oP- The "correction" 
operator then behaves as 

C - O (Vn) , (21) 

and we have no true control over the corrections in the 
mean-field expansion which ought to be negligible. The 
mean-field approximation can only be consistent (if at 
all) with a modified GP equation (|15|l inducing a renor- 
malization of the coupling according to Eq. (|16|) . How- 
ever, even given these modifications, the applicability 
of the mean-field expansion is not obvious since higher- 
order operator products can yield 0{N) contributions af- 
ter the fc-summation/integration and thus the hierarchy 
C ^ X is not evident. Similarly, the first-order correction 
in Eq. p6ll is comparable to the zcroth order (i.e., of the 
same order in N), which hints at the fact that all orders 
must be taken into account in a suitable way. 

It is illuminating to compare the employed dilute-gas 
limit {N 1 oo particles in a fixed volume with gN remain- 
ing constant) with the thermodynamic limit (N "f oo par- 
ticles in an increasing volume V with g and N/V remain- 
ing constant): Translation of the scaling in Eq. H12|l to 
the thermodynamic limit yields an interaction potential 
whose range increases proportional to the system size V 
whereas its strength decreases accordingly. It is not very 
surprising that the mean-field expansion is very good in 
this case. The analogue of Eq. H18|l in the thermodynamic 
limit is an interaction potential with constant strength 
and range (where the applicability of the mean-field ex- 
pansion is less obvious). Finally, Eq. H2U|) corresponds to 
a potential with decreasing range and increasing strength 
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(in the thermodynamic hmit) . Again, is not very surpris- 
ing that such a scahng might generate difficulties in the 
ordinary mean-field expansion and requires taking into 
accout all orders in a suitable way. 



an instability, cf. [17| ■ Further expectation values such as 
the total energy E = (H) can be calculated analogously 

V. CONCLUSION 



IV. DIPOLE-DIPOLE INTERACTION 

As an example for interactions with a finite range, let 
us consider a dipole-dipole force cx gd in addition to the 
usual contact repulsion oc g. Dipole-dipole interactions 
between atoms [l^ can either be induced by an external 
electric field or be due to an intrinsic magnetic dipole 
moment. E.g., for magnetic dipoles, gd — /xgd^/S; Bose- 
Einstein condensation of chromium, which has a ground 
state moment of dm = 6/iB, has been achieved recently 

If the dipole moments of all atoms/molecules are 
aligned along the z-axis, the dipole-dipole interaction po- 
tential reads (in three spatial dimensions) 



Vddir) 



3gd 1 - 3zV|r| 
An IrP 



(22) 



In addition to this long-range interaction, the particles 
are subject to a short-range repulsion whose impact can 
be represented by a contact potential cx S^{r — r'). Con- 
sequently, the Fourier transformed potential reads for in- 
termediate momenta 



Vir^tik) = g + gd [ -ri - I 



3_fc? 

fc2 



(23) 



The (ideal) dipole-dipole interaction potential behaves 
as which is not integrable and, strictly speaking, is 
therefore just at the limit of applicability of our analysis - 
for example, Vint is not well-defined at fc = 0, which com- 
plicates the introduction of an effective coupling g (finite- 
size effects etc.) Nevertheless, we may regard the non- 
integrable dipole-dipole interaction potential as a limit 
of integrable potentials and calculate the corresponding 
corrections beyond the s-wave scattering approximation. 
For example, the quantum depletion is modified from the 
pure contact case via 



3/2 

247r2 



\/%d + g (gd + 5g) 




This expression is valid for g > gd only, since otherwise 
the excitation spectrum (for a homogeneous condensate) 
and hence also (x^x) contain imaginary parts indicating 



Based on the point-particle approximation of Eq. Q , 
we derived the mean-field expansion Eq. Q for dilute 
Bose-Einstein condensates with arbitrary weak finite- 
range two-particle interaction potentials, obeying suit- 
able scaling behavior in the large- limit. It turns out 
that, although the gas rapidly becomes ever more dilute 
in the large A^ limit (the gas parameter (g'^g)^^^ oc 1/A^), 
the validity of the mean-field approximation strongly de- 
pends on the detailed scaling behavior of the particle 
interaction potential. Therefore, care needs to be exer- 
cised in applying the mean-field approximation ~ if one 
does not take the detailed structure of the bare interac- 
tion potential into account, one possibly encounters in- 
consistencies with the basic assumptions the mean-field 
expansion is built on. 

Apart from exploring limits of the mean-field approx- 
imation, our derivation facilitates the calculation of the 
impact of microscopic details of the interaction beyond 
the s-wave scattering approximation ITij . Another ad- 
vantage of the presented approach is the natural emer- 
gence and ab initio derivation of the cut-off fccut as a mi- 
croscopic property of the interaction potential - instead 
of a cut-off fccut introduced ad hoc for the regularization of 
pseudo-potentials cx S'^{i — r'), see, e.g., j20|. However, 
one should bear in mind that the presented method re- 
quires weak potentials - if nonperturbative effects such as 
bound states or total reflection at a finite radius become 
important, the direct mean-field ansatz ^ cannot be ap- 
plied in this way. Instead of homogeneous plane waves, 
one has to start from atomic/molecular eigenfunctions 
in this case and the ansatz (|3Jl can only work as a low- 
energy effective description. Nevertheless, if it is possi- 
ble to divide a nonperturbative Vint{r — r') into a strong 
ultra-short-range part and a comparably weak remain- 
ing contribution which is more spread out, the presented 
analysis should be applicable to the latter part, provided 
the ultra-short-range part can be replaced in an adequate 
manner by an integrable (pseudo-)potential. An exam- 
ple for the application of our analysis to the longer-range 
contribution we gave in the last Section, using a dipole- 
dipole interaction potential. 
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